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ABSTRACT 

We report the discovery of an extrasolar planet, KMT-2015-lb, that was detected using the microlensing 
technique. The planetary lensing event was observed by KMTNet survey that has commenced in 2015. With 
dense coverage by using network of globally distributed telescopes equipped with very wide-held cameras, the 
short planetary signal is clearly detected and precisely characterized. We hnd that KMT-2015-lb is a giant 
planet orbiting a low-mass M-dwarf host star. The planet has a mass about twice that of Jupiter and it is located 
beyond the snow line of the host star. With the improvement of existing surveys and the advent of new surveys, 
future microlensing planet samples will include planets not only in greatly increased number but also in a wide 
spectrum of hosts and planets, helping us to have a better and comprehensive understanding about the formation 
and evolution of planets. 

Subject headings: gravitational lensing; micro - planetary systems 


1. INTRODUCTION 

Since t he hrst discovery b y [Wols zczan & Fraill (119921) fol¬ 
lowed by iMavor & OuelozI (1 19951) . many exo-planets have 
been discovered. With the Kepler mission, the number of 
known planets explosively increased and now reaches ^2000. 
Most of them were discovered by either the transit or radial- 
velocity methods. 

Although relatively few, planets have also been discov¬ 
ered using the microlensing method. Due to the fact that 
these planetary systems are detected through their gravita¬ 
tional fields rather than their radiation, this method makes 
it possible to detect planets around faint stars and even dark 
objects. Furthermore, microlensing is sensitive to planets on 
wide orbits beyond the snow line, which separates regions of 
rocky planet formation from regions of icy planet formation, 
while other major planet detection techniques are sensitive to 
close-in planets. Being able to detect planets that are difficult 
to be detected by other techniques, the method is import ant for 
the co mprehensive understanding of planet formation (iGaudil 

One of the important reasons for the small number of mi¬ 
crolensing planets is the difficulty of observation. The plan¬ 
etary signal of a microlensing planet is a short-term pertur- 
bati on to the smooth lensing light curve induced b y the host 
star (iMao & Paczvnskill99itlGould & Loebll 19921) . To detect 
short planetary signals, the previous generation of lensing ex¬ 
periments applied a strategy where lensing events were de¬ 
tected by wide-held surveys and events detected by surveys 
were in tensively monitored using multiple narrow-held tele¬ 
scopes (lAlbrow et al.lll998h . However, the limited number of 
follow-up telescopes makes it difficult to intensively monitor 
all alerted events, which number reaches several hundreds at 
a given moment. As a result, the detection efficiency of mi¬ 


crolensing planets has been low. 

In this paper, we report the discovery of a giant planet 
KMT-2015-lb orbiting a low-mass M dwarf. The planetary 
lensing event was observed by a new survey that achieves 
round-the-clock coverage of lensing events with a high ca¬ 
dence by using global network of telescopes equipped with 
very wide-held cameras. 

The paper is organized as follows. In Section 2, we de¬ 
scribe the observation of the planetary microlensing event by 
the new lensing survey, including instrument, data acquisition 
and reduction process. In Section 3, we give a description 
of the modeling procedure conducted to analyze the observed 
lensing light curve. We provide the estimated physical param¬ 
eters of the discovered planetary system in Section 4. Finally, 
we summarize the result and make a brief discussion about 
the result in Section 5. 

2. OBSERVATION 

The planet was discovered from the observation of a 
microlensing event that occurred on a star located to¬ 
ward the Galactic bulge held. The equatorial coor¬ 
dinates of the lensed star (source) are (RA,DEC) = 
(17*'58”39**.01,-28°01'54".l), that correspond to the Galac¬ 
tic coordinates {l,b) = (2.24°,-2.00°). The lensing-induced 
brightening of the s ource star was notic ed in early February, 
2015 by the OGLE (lUdalski et al.ll2015h survey and lasted for 
about a month, during which the bulge held can be seen for 
less than an hour per night. 

The event was also observed by the KMTNet (Korea Mi¬ 
crolensing Telescope Network) lensing survey that started its 
test observation in February, 2015, which matches the occur¬ 
rence time of the event. The event was dubbed as KMT-2015- 
BLG-0048 in the KMT event list. The survey uses three iden¬ 
tical telescopes that are located at Cerro Tololo Interameri- 
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Fig. 1.— Light curve of the microlensing event KMT-2015-BLG-0048. Solid and dashed curves represent the best-fit models from binary (planetary) and 
single-mass lens modeling, respectively. 


can Observatory, Chile (KMT CTIO), South African Astro¬ 
nomical Observatory, South Africa (KMT SAAO), and Siding 
Spring Observatory, Australia (KMT SSO). At the time of the 
event, KMT SSO was not online and the event was observed 
by two telescopes, KMT CTIO and KMT SAAO. Each tele¬ 
scope has a 1.6m aperture and is equipped with a mosaic cam¬ 
era composed of four 9Kx9K CCDs. Each CCD has a pixel 
size of 10 microns corresponding to 0.4 arcsec/pixel and thus 
the camera has 4 deg^ field of view (iKim et al.ll2015h . With 
the global network of wide-held telescopes that are separated 
by ^8 hrs in longitudes, the survey continuously observes the 
Galactic bulge held with a cadence of 1 observation every 
10 minutes. Images were taken in Johnson-Cousin I and V 
bands. Erom the observation of the event, we acquired 786 
/-band and 54 V-band images from KMT CTIO and 1117/- 
band images from KMT SAAO observations. We note that 
/-band data are used to analyze the lensing event and V-band 
data are used to constrain the source star. 

Photometry of the images were conducted using a cus¬ 
tomized pip eline that is b a sed on the difference image anal¬ 
ysis code of lAlbrow et al.l (l2009l) Since data are taken by the 
different telescopes and photometry based on an automated 
pipeline often results in inaccurate error estimation, we renor¬ 
malize error bars of the individual data sets by 

o•' = k(c^g-l-cr^J‘/^ (1) 

where ctq is the error estimated from the pipeline, is a 
factor used to make the cumulative distribution function of 
as a function of lensing magnihcation linear, and k is a scaling 
factor to make x^ per degree of freedom (dof) become unity. 
We note that the factor Cniin is required to ensure that each 
data set is fairly weighted according to its error bars. 

In Eigure[T] we present the light curve of KMT-2015-BLG- 
0048. Compared to the continuous and symmetric light curve 
of a single-mass event, the light curve exhibits a short-term 
perturbation during 7080.0 < HJD-2450000 < 7082.5. The 


TABLE 1 

Lensing parameters 


Parameters 

Values 

X^/dof 

1896/1896 

to (HID) 

2457083.081 ±0.004 

WO 

0.221 ± 0.003 

tE (days) 

10.90 ±0.10 

s 

0.955 ± 0.003 

g(io-3) 

6.84 ±0.18 

a (rad) 

5.352 ± 0.003 

p(10-3) 

44.7 ± 0.8 


perturbation shows a feature that is composed of a depression 
centered at HJD-2450000 ^ 7081.5 and brief bumps at both 
edges of the depression. Such dips, usually surrounded by two 
bumps, are a generic feature of lensing systems with small 
mass ratios q <^1 with normalized planet-star separations s < 
1, i.e., planets inside the Einstein ring. This is because the 
star’s gravity generates two images, one inside and the other 
outside the Einstein ring. The former, being a saddle point on 
the time delay surface, is easily suppressed if the planet lies 
in or near the path, thereby causing relative demagnification, 
and hence a dip in the light curve. To be noted is that the 
major structure of the anomaly feature was well covered by 
the survey observation despite the short time window each 
night toward the field. 

3. ANALYSIS 

Keeping in mind that the anomaly pattern is likely to be 
produced by a binary lens with a low mass ratio, we conduct 
binary-lens modeling. Basic binary-lens modeling requires at 
least 9 parameters, including 7 principal parameters for the 
lensing system and 2 flux parameters for each observatory. 
The first 3 of the principal parameters describe the source ap¬ 
proach with respect the lens, including the time of the closest 
source approach to the reference position of the lens, to, the 
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Fig. 2.— Lens system geometry. The upper panel shows the source trajec¬ 
tory (straight line with an arrow) with respect to the lens components (marked 
by M{ and M 2 ) and caustics (closed concave curve) and the lower panel shows 
the light variation with the progress of the source position. Lengths are scaled 
to the Einstein radius and the source trajectory is aligned so that the progress 
of the source matches the light curve shown in the lower panel. The inset in 
the upper panel shows the wide view and the major panel shows the enlarged 
view around the caustic. The empty circles on the source trajectory repre¬ 
sent the source positions at the times of observation and the size indicates the 
source size. The dotted curve in the lower panel is the light curve expected 
for a point source. 


lens-source separation at that moment, mq (impact parameter), 
and the time for the source to cross the angular Einstein ra¬ 
dius 0E of the lens, (Einstein time scale). Eor the reference 
position of the lens, we use the center of mass. The other 
3 principal parameters describe the lens binarity including the 
projected separation s and the mass ratio q between the binary 
components, and the angle between the source trajectory and 
the binary axis, a. We note that the impact parameter mq and 
the binary separation s are normalized by 9^. The other pa¬ 
rameter defined as the ratio of the angular source radius to the 
Einstein radius, p = is needed to describe light curve 

deviations affected by finite-source effects. 

Eor some lensing events, the basic parameters are not ade¬ 
quate to precisely describe observed lensing light curves. The 
known cause s of such deviations includ e the parallax effect 
jGould||1992h and the lens orbital effect (lAlbrow et al.ll2000t 
lAn et ai.l2Q02HJung et al.l2013h . The parallax effect is caused 
by the deviation of the observer’s motion from rectilinear due 
to the orbital motion of the Earth around the Sun. On the other 
hand, the lens-orbital effect is caused by the orbital motion of 
the lens. Such effects are important for long time-scale events 
for which the orbital motion of the Earth and/or lenses are 
important. Eor KMT-2015-BLG-0048, which lasted for ^1 
month, we find that these higher-order effects are negligible. 

Modeling the light curve proceeded in several steps. Eirst, 
we conduct thorough a grid search for solutions in the param¬ 
eter space of (s, q, a), for which lensing light curves vary sen¬ 
sitively to the change of the parameters. In this process, other 



Fig. 3.— Source position in the instrumental color-magnitude diagram of 
nearby stars with respect to the centroid of the giant clump. 


parameters are searched for by using a downhill approach. 
Second, we investigate all possible local solutions in the pa¬ 
rameter space in order to check the existence of degenerate 
solutions where different combinations of the lensing param¬ 
eters result in a similar light curve. Einally, we search for the 
global solution by comparing values of the identihed lo¬ 
cal solutions. Eor the downhill minimization, we use the 
Markov Chain Monte Carlo (MCMC) method. 

Erom the search for a solution, we hnd that the event was 
produced by a planetary system where the planet/host mass 
ratio is q = 6.84 x 10“^ and the projected planet-host sepa¬ 
ration, s = 0.96, is slightly less than the Einstein radius of 
the lens. In Table [T] we present the best-ht lensing param¬ 
eters, where the error bars are estimated from the scatter of 
the MCMC chain. In Eigure |2] we present the geometry of 
the lens system corresponding to the solution. Due to the 
resonance of the projected separation to 9e, i.e. i ^ 1, the 
lens system forms a single large caustic around the host of 
the planet. The source passed the backside of the atTowhead- 
shaped caustic. The depression in the light curve occurred 
when the source was in the demagnification valley between 
the two protrudent cusps that caused the brief bumps on both 
sides of the depression. The source crossed the tip of one of 
the cusps during which the light curve shows a clear finite- 
source signature from which we accurately measure the nor¬ 
malized source radius p. To be noted is that the measured 
value of p = 44.7 x 10“^ is significantly larger than typical 
values ^10“^ for main-sequence and ^10“^ for giant source 
stars. Eor a given size of a bulge star, this suggests that the 
Einstein radius is very small. 

Since finite-source effects are clearly detected, it is possi¬ 
ble to determine the angular Einstein radius. The Einstein 
radius is determined by 0^ = 9^!p. The normalized source ra¬ 
dius p is measured from the analysis of the light curve around 
the planetary perturbation. The angular radius of the source 
star, 0*, is determined from its de-reddened color (y-/)o 
and brightness /q. To measure color and brightness, we lo¬ 
cate the source star in the instrumental(uncalibrated) color- 
magnitude diagram of neighboring stars in the same held and 
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calibrate the color and brightness based on the location of the 
centroid of the giant clump (GC), for which its de-reddened 
color (V-/)o,GC and brightness /q gc are known to be con¬ 
stant and thus can be used as a standard candle. Figure [2 
shows the locations of the source and giant clump centroid 
in the instrumental color-magnitude d iagram. By adopting 
{V - /)o,GC = 1 06 (iBensbv et al.l 1201 Ih and brig htness /n or 
accou nting for variation with Galactic longitude (iNataf et al.l 
1201 3h . we find that (V—I,I)o = (0.78,14.5), indicating that the 
source is a G-type giant star. We then convert V- I into V -K 
using the color-color relation of iBessell & Bred (119881) and 
finally determine 0^ u sing the color-angular radius relation of 
iKervella et al.l ( 120041) . We find that the angular source radius 
is 0* = 4.39 ± 0.56 /ras. Then the angular Einstein radius is 

6»e = 0.098 ±0.013 mas. (2) 

Combined with the Einstein time scale estimated from lens 
modeling, the relative lens-source proper motion is deter¬ 
mined as 

^ =3.29±0.43masyr-'. (3) 

Ie 

As expected, the estimated Einstein radius is substantially 
smaller than ^ 0.5 mas for typical Galactic microlensing 
events. 

4. PHYSICAL PARAMETERS 

Eor the unique determination of the mass M and distance 
Di to the lens, it is required to simulta neously meas ure both 
the lens parallax and the Einstein radius (iGouldl 19921) . Never¬ 
theless, the above measurements of the Einstein radius 9e and 
the proper motion ^ can be combined with a Galactic model to 
statisti cally constrain M and Di yia a Bayesian analysis. We 
use the lHan & Gouldl ( |1995h model, whose matter distribution 
is based on a double-exponential disk and a triaxial bulge. The 
disk yelocity distribution is assumed to be Gaussian about the 
rotation yelocity and the bulge yelocity distribution is a tri¬ 
axial Gaussian with components deduced from the flattening 
yia the tensor yirial theorem. The mass function is based on 
Hubble Space Telescope star counts. 

The results of this Bayesian analysis 

M = O.18±O.12M0; = 8.2± 0.9 kpc (4) 

can be easily understood based on simple physical reasoning. 
Eirst, the Enstein radius is by definition the root-mean-square 
to the lens mass and the lens-source relatiye parallax TTrei, i.e. 


6'e= V/^MTTrei; ^ = ^TTT - 1 ■ (5) 

c^AU Mq 

Because 9e = 0.098 mas is unusually small, either M or TTrei 
should be small: (M/MQ)(7rrei//tas) = 1.2. Thus, if the lens 
is a star (M > 0.08 Mq), it is Tr^ei < 15 pas and so lies in the 
Galactic bulge. This, combined with the fact that the Galac¬ 
tic density profile along the line of sight peaks strongly in 
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the bulge, heayily fayors bulge lenses. Einally, the measured 
proper motion is yery typical of bulge lenses, but only half the 
size of typical proper motions. 

In principle, the lens could be a star of any mass in the 
bulge. Howeyer, the measurement of 9e telles us that if it 
were a solar mass star, then TTrei =1.2 pas, which would imply 
a distance from lens to source of only Ds-Di = 70 pc. It is 
this small amount of ayailable phase space for heayier lenses 
that driye the Bayesian estimate toward low mass lenses. 

The measurement of q from the lens model then directly 
yields an estimate of the planet mass 

Mp = qM = 2.2±lAMj, (6) 

while the measurement of s yields the projected lens-host sep¬ 
aration 

£/j_=iDz,6»E = 0.76 ± 0.08 AU. (7) 

Considering the low temperature of the host star, the planet is 
located beyond the ice line. 

5. SUMMARY AND DISCUSSION 

We reported the discoyery of an extrasolar planet KMT- 
2015-lb that was detected by using the microlensing tech¬ 
nique. The planetary lensing eyent was obseryed by the KMT- 
Net suryey that has commenced in 2015. Despite the short 
time window in the early bulge season, unambiguous detec¬ 
tion and precise characterization of the planetary system was 
possible due to the dense coyerage of the planet-induced per¬ 
turbation by using network of globally distributed telescopes 
equipped with yery wide-field cameras. We find that KMT- 
2015-lb is a giant planet orbiting a low-mass M-dwarf host 
star. The planet has a mass about twice that of the Jupiter and 
it is located beyond the snow line of the host star. 

Cool M dwarfs far outnumber sun-like stars and thus under¬ 
standing the process of planet formation around them is im¬ 
portant. Eurthermore, small masses and low luminosities of 
M dwarfs proyide leyerage on conditions of planet formation, 
enabling to check the yalidity of existing formation theories 
and refine suryiyed theories. With the improyement of exist¬ 
ing suryeys and the adyent of new suryeys, future microlens¬ 
ing planet sample will include planets not only in greatly in¬ 
creased number but also in wide spectrum of hosts and plan¬ 
ets, helping us to haye better and comprehensiye understand¬ 
ing about the formation and eyolution of planets. 
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